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ABSTRACT 

Development of satisfactory techniques for detecting 
change in coastal zone environments is required before 
operational monitoring procedures can be established. In 
an effort to meet this need a study was directed toward 
developing and evaluating different types of change detec- 
tion techniques, based upon computer-aided analysis of 
Landsat multispectral scanner (MSS) data, to monitor these 
environments . 

The Matagorda Bay estuarine system along the Texas 
coast was selected as the study area. Land use within 
this area is divided principally among agriculture, range- 
land, urban, industry, recreation and large marsh-covered 
tracts . 

Four change detection techniques were designed and 
implemented for evaluation: a) post classification com- 

parison change detection, b) delta data change detection, 
c) spectral/temporal change classification, and d) layered 
spectral/temporal change classification. Each of the four 
techniques was used to analyze a Landsat MSS temporal data 
set to detect areas of change of the Matagorda Bay region. 

The post classification comparison technique reliably 
identified areas of change and was used as the standard for 
qualitatively evaluating the other three techniques. The 
layered spectral/temporal change classification and the 
delta data change detection results generally agreed with 
the post classification comparison technique results; how- 
ever, many small areas of change were not identified. 

Major discrepancies existed between the post classification 
comparison and spectral/temporal change detection results. 


INTRODUCTION 

Our nation's coastal zones are locales of unending change. Effective in- 
ventorying and monitoring of these changes are required for resource management. 
Conventionally, monitoring of these areas has been accomplished either by local 
on-site surveys or photointerpretation of aerial photographs, and change detec- 
tion has required the manual comparison of various single date analyses. For 
large area inventories, both local surveys and photointerpretative techniques 
require large data collection efforts and are time consuming and subjective 
in nature. 


♦This work was accomplished under the National Aeronautics and Space Adminis- 
tration Contract No. NAS9-14016. Journal Paper No. 6689 Purdue University, 
Agricultural Experiment Station. 
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Consequently, it would be desirable to develop and implement a computer- 
ized change detection procedure suitable for large area inventories. (1) How- 
ever, before operational monitoring procedures can be established, satisfactory 
computer-aided change detection techniques must be developed. In order to meet 
this need, this study was initiated toward developing and evaluating various 
change detection techniques based upon computer-aided analysis of Landsat multi- 
spectral scanner (MSS) data to monitor coastal zone environments. 

STUDY SITE 

A portion of the Matagorda Bay estuarine system located along the Texas 
coast was selected for this study. This area is characterized by a great di- 
versity in geography, resources, climate, waterways and estuaries similar to 
the entire coastal zone of Texas. Also, it is subjected to numerous natural 
hazards such as shoreline erosion, land surface subsidence, stream flooding, 
hurricane tidal surges and active surface faulting. Lower elevations of the 
coastal area are covered with marshes and swamps, and the standplains are 
chiefly sandy materials. The coastal upland areas consist mainly of clay and 
mud materials deposited as sediments during the Pleistocene era. Land use 
within this study site is divided principally among agriculture, rangeland, 
urban, industry, recreation and large marsh-covered tracts. (2) 

APPROACH 

Data 


Landsat MSS data collected on November 27, 1972 and February 25, 1975, 
were used as the principal data sources for this study. These data were geo- 
metrically corrected (i.e., rotated, deskewed and rescaled) (3) , overlaid and 
precision registered to ground control points selected from appropriate U.S. 
Geological Survey (USGS) 1 \ minute topographic quadrangle maps. This proce- 
dure produced a multidate eight-channel data set at a scale of 1:24,000 and 
registered points in the data to their exact ground position. Reference data 
utilized to support the analysis of the Landsat data included 1970, 1971 and 
1975 color and color infrared aerial photography. Geologic Atlases compiled by 
the Texas Bureau of Economic Geology, and Spectral Environmental Classification 
overlays produced by Lockheed Electronics Corporation. 

Single Date Analysis 

A non-supervised approach utilizing the maximum likelihood classifier as 
implemented by LARSYS (4) was used to produce a classification for each date 
for areas corresponding to individual USGS minute quadrangle maps. The 
resultant classifications included informational categories of agriculture, 
rangeland, grasslands, wooded areas, swamp timber, fresh and brackish marshes, 
mud and sand flats, beach and beach ridges, residential areas, and open water 
classes. This initial study demonstrated that digital analysis of Landsat MSS 
data could be used to inventory the coastal environment effectively. (5) 

Change Detection Analysis 

Four specific test sites were selected within the Matagorda Bay estuarine 
system for implementation and evaluation of change detection techniques. These 
were the areas included in the Austwell, Port Lavaca-E, Port O'Connor and Pass 
Cavallo-SW USGS minute topographic quadrangles. 

Four techniques for detecting change were designed and implemented for 
evaluation (6) : a) post classification comparison, b) delta data change detec- 

tion, c) spectral/temporal change classification, and d) layered spectral/ 
temporal change. All four techniques require registration of the data prior 
to analysis for change. A brief description of each technique follows. 
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Post Classification Comparison (7) 


Post classification change detection is based on the comparison of inde- 
pendently produced spectral classifications. The comparator is simply a pro- 
cessor that "compares" two standard LARSYS classifications utilizing class 
pairs specified by the analyst and generates a results tape indicating areas 
of change. Thus, the change classes derived are defined by the analysts 
rather than being identified spectrally using pattern recognition techniques. 
By properly coding the classification results for times ti and t 2 , the analyst 
can produce change maps which show a complete matrix of changes. In addition, 
selective grouping of classification results allows the analyst to observe any 
subset of changes which may be of interest. 


Delta Data Change Detection (8) 


The delta data change detection technique is based upon the classification 
of a multispectral difference data set. A delta (subtraction) transformation 
combines two n-channel multispectral data sets obtained at different times and 
produces a multispectral delta data set having n-channels. Since the LARSYS 
system assumes all data samples used are non-negative, a bias is added to each 
difference so that the resulting delta data is non-negative. Thus, the com- 
plete transformation is simply: 


Where : 



ti 

t 2 


Multispectral value for channel k; i,j = 1 

N x N image 

Row 

Column 

Channel 

Bias for Channel k 

First date 
Second date 


.N assuming an 


This approach assumes that specific changes from one time to another will pro- 
duce a non-negative result which can be detected either by examining the 
images directly or by first classifying the delta data, then examining the 
results displayed in inage form. The classification method is similar to that 
for classification of ordinary multispectral data; however, it emphasizes clas- 
sification of multispectral change instead of changes in multispectral classi- 
fication. 


Spectral/Temporal Change Classification 

This technique is based upon a single analysis of a multidate data set to 
identify change. In a two-date Landsat data set, eight channels of data would 
be analyzed using standard pattern recognition techniques. By using data sets 
collected under similar conditions at nearly the same day of the year but from 
different years, one would expect that classes where change is occurring would 
have statistics significantly different from where change had not occurred and 
could be so identified. 

Layered Spectral/Temporal Approach (9,10,11) 

The layered or decision tree classifier is essentially a maximum likeli- 
hood classifier using multi-stage decision logic. It is characterized by the 
fact that an unknown sample can be classified into a given class using one or 
more decision functions and channels in a successive manner. This classifi- 
cation strategy can be most easily illustrated by a tree diagram consisting 
of a root node, a number of non-terminal nodes or decision stages (layers) and 
terminal nodes (i.e., the decision making procedure terminates with the 
unknown sample being assigned to the class at a terminal) . A non-terminal 
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node is an intermediate decision, its immediate descendent nodes representing 
the possible outcomes of that decision. 

To specify a decision tree uniquely, two sets of information are neces- 
sary: a) how the terminal and non-terminal nodes are linked and b) the deci- 

sion functions and channels of all the non-terminal nodes. The decision tree 
can be constructed manually or by an automatic optimized logic tree design 
procedure. For the purpose of change detection a hybrid of these two proce- 
dures was utilized. Decision trees for each of the two dates, ti and t 2 » were 
obtained automatically and then manually linked introducing within the tree a 
logic for detecting the desired changes. 

RESULTS 

A quantitative evaluation of the change detection procedures developed 
would require the collection of a comprehensive set of reference data (ground 
observations and/or aerial photography) coincident with the Landsat overpass. 
However, coincident reference data were not available for this particular pro- 
ject (i.e., reference data for the 1972 Landsat data consisted of aerial pho- 
tography collected in 1970, 1971). Thus, only a qualitative evaluation of the 
results of the change detection techniques could be accomplished. 

Since the comparison of results from two single date classifications is 
the basis of the post classification comparison change detection procedure, 
the results of this technique were used as the standard for evaluating the 
results from the other three procedures. Correlation of the unsupervised spec- 
tral classes with informational classes was accomplished by: 

a) associating the spectral classifications with aerial photography 

when appropriate, v 

b) utilizing a ratio, A = y r calculated for each spectral class, where V 

is the relative intensity of the mean spectral values of the 

visible wavelengths {(0.5 to 0.6ym) + (0.6 to 0.7ym)} and IR is the 
relative intensity of the mean spectral values of the reflective 
infrared wavelengths {(0.7 to 0.8um) + (0.8 to l.lym)} and 

c) by summing the relative intensity values ("summed response") of all 
four bands to determine the magnitude of relative spectral responses 
for each spectral class. 

By observing the aerial photography, the ratio A and the summed response, the 
analyst delineated major vegetation and land use categories within the coastal 
zone area. Results from the Port O'Connor quadrangle, which are representative 
of the four test areas, are presented for illustrative purposes. 

Utilizing the post classification comparison technique, the November 1972 
and February 1975 dates for the Port O'Connor test site were classified into 
25 and 23 spectral classes, respectively. These classes were grouped into 6 
ecological categories: a) urban, b) woody/herbaceous vegetation, c) submerged 

vegetation and tidal flats, d) spoil areas, e) water and f) burned areas. The 
logic comparator was constructed to identify the class pairs listed in Table I. 
The map produced by this comparison (Figure 1) illustrates the areas of change. 

Cluster analysis of the delta data set produced 13 delta spectral classes 
representing seven identifiable informational classes (Table II) . These infor- 
mation classes were not in a 1:1 correspondence with the post classification 
comparison classes but rather were dispersed over a number of the classes. 
However, visual inspection of the delta data change results map showed a good 
general agreement between the post classification comparison and delta data 
change results in that the no-change water/land interfaces and the delineation 
of burn areas were in close agreement. Pixels classified into the mixed and 
urban change classes were widely dispersed through the land and marsh areas. 
This situation may be related to the inherent difficulty of using present 
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Landsat resolution and wavelength bands to discriminate reliably among certain 
surface features and not to the change detection procedure. 

To evaluate the spectral/temporal approach to detecting change, three 
areas from the eight-channel data set were clustered to identify potential 
change classes. These included the agricultural region, intercoastal waterway 
and the barrier island, areas where change was expected to be occurring. Sta- 
tistics developed during clustering were used to identify potential change 
classes using the following empirically-derived criteria: 

a change class was probable if the ratio (A = for each of the two 

dates differed by more than 0.30, and the summed response differed 

by more than 25.0. 

Fifty-two spectral classes were identified resulting in five information 
classes (Table III) , four of which represented change. 

Comparison of the results of the spectral/temporal approach to the post 
classification comparison technique indicated major disagreement. Pixels clas- 
sified into the woody to woody class by the post classification comparison 
technique were classified as vegetation to soil, and soil to vegetation change 
classes (Table III) by the spectral/temporal change detection technique. Sim- 
ilarly, the pixels classified into the submerged-submerged, any non-submerged- 
submerged and water-water classes were classified as a water to soil change 
class. Only the burn change class agreed well with the post classification 
comparison change results. 

Detection of change using the layered spectral/temporal approach involved 
developing a layered decision tree based upon time: statistics for time: (No- 
vember 1972) and attaching a layer decision tree based upon time 2 statistics 
for time* (February 1975) to the appropriate terminal nodes of time: producing 
a "change decision tree". Results of this procedure produced 47 spectral 
classes, and 18 informational classes, 10 of which were change classes (Table 
IV) . 


Results of the layered spectral/temporal approach were in good agreement 
with those obtained by the post classification comparison procedure. Burn 
areas and the no-change land/water interfaces were very similar. Urban areas 
were well identified; however, occasionally the dredged spoil areas along the 
intracoastal waterway were placed into the urban class due to the similarities 
in spectral response. Some of the confusion between classes was attributed to 
the trimming of the decision trees (which was necessary to accomodate the tree 
in computer memory) . 


CONCLUSIONS 

In devising a change detection procedure one property to consider is the 
relative complexity of the method; all other things being equal, a simpler 
approach would be more desirable for implementation reasons. However, more 
complex methods may have more performance potential in the long run. On a 
scale of increasing complexity the four methods investigated may be ordered 
as follows: 

1. Delta data method . This method requires only a simple subtraction 
followed by a single classification. 

2. Post classification comparison . This method requires two separate 
classifications followed by a logical comparison. 

3. Spectral/temporal classification . While this method requires only a 
single classification, it is a vastly more complex one, requiring 
more classes and probably more features. 

4. Layered spectral/temporal classification . This method involves not 
only a complex classification but also a priori knowledge of the 
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logical interrelationship of the classes as well. 


In a change detection analysis, change detection error can arise from 
classification error at date one, date two or both times. To evaluate pre- 
cisely the results of various change detection procedures, comprehensive 
ground information must be available to identify any error condition. 

Since adequate reference data for a thorough evaluation was not available, 
the post classification comparison change detection results were used as the 
standard for qualitative evaluation of the results from the other three proce- 
dures. Areas of major change were reliably identified by the post classifica- 
tion comparison technique. This is due in part to the fact that the analysis 
procedures required are routine and quite well understood. 

Evaluation of the other three change detection techniques indicates the 
following: 

a) The delta data change detection method may be too simple to deal 
adequately with all the factors involved in detecting change in a natural 
scene. Too much information may be discarded from the data in the subtraction 
process whereby only the four band difference data remains from the two sets 
of original four bands. Images created from the delta data may be quite 
useful, however, in qualitatively assessing change by image interpretation. 

b) The spectral/temporal and layered spectral/temporal change detection 
methods cannot be ruled out at this point as each appears to have undeveloped 
potential. The layered spectral/temporal change detection results in par- 
ticular showed best agreement with the post classification comparison results. 
More complex methods usually require more carefully drawn data inputs together 
with greater user understanding to achieve their potential. 
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TABLE I. CLASS PAIRS DEFINED FOR THE POST CLASSIFICATION 
COMPARISON CHANGE DETECTION TECHNIQUE FOR TWO 
DATA SETS FOR THE PORT O'CONNOR QUADRANGLE. 


Class No. 


Informational Classes 
NOV 72 FEB 75 


1 

2 

3 

4 

5 

* 6 

* 7 

* 8 
* 9 
*10 
*11 


urban 

woody/herbaceous 
spoil 
submerged 
water 
any non-burn 
any non -spoil 
any non-woody 
any non-water 
any non- submerged 
any non-urban- 


* change classes 


-urban 

-woody/herbaceous 

-spoil 

-submerged 

-water 

-burn 

-spoil 

-woody 

-water 

-submerged 

-urban 


TABLE II. CLASSES DELINEATED BY DELTA DATA CHANGE DETECTION 
METHOD AND THEIR CORRESPONDING POST CLASSIFICATION 
COMPARISON INFORMATIONAL CLASSES FOR THE PORT 
O'CONNOR QUADRANGLE. 

Class No. Informational Classes 


* 1 
* 2 

* 3 

* 4 

* 5 


6 

7 

8 
9 

10 

11 

*12 

13 


any non-spoil to spoil 
any non-spoil to spoil 

mixed and any non-submerged to submerged 
mixed and any non-urban to urban 
woody to wood and any non-submerged to 
sumberged and any non-urban to urban 
woody to woody 

woody to woody and water to water 

woody to woody 

water to water 

water to water 

water to water 

any non-burn to burn 

water to water 


* change classes 


TABLE III. CLASSES DELINEATED BY SPECTRAL/TEMPORAL 

CHANGE DETECTION METHOD FOR PORT O'CONNOR 
QUADRANGLE . 


Class No. 


Informational Classes 


*1 

*2 

*3 

4 

*5 

6 


Vegetation to soil 
Soil to vegetation 
Vegetation to burn 
No change land class 
Water to soil 
No change water class 


* change class 
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TABLE IV. CLASSES DELINEATED BY LAYERED SPECTRAL/TEMPORAL 
CHANGE DETECTION METHOD FOR THE PORT O'CONNOR 
QUADRANGLE . 

Class No. Informational Classes 



NOV 72 

FEB 75 

* 1 

woody 

urban 

2 

urban 

urban 

3 

woody 

woody 

* 4 

urban 

woody 

* 5 

woody 

burn 

* 6 

urban 

burn 

7 

water 

water 

* 8 

spoil 

water 

* 9 

submerged 

water 

*10 

water 

spoil 

11 

spoil 

spoil 

*12 

submerged 

spoil 

*13 

water 

submerged 

*14 

spoil 

submerged 

15 

submerged 

submerged 

16 

water 

confusion 

17 

spoil 

confusion 

18 

submerged 

confusion 


* change classes 
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